Spin thermal conductivity of the Haldane chain compound Y 2 BaNiC>5 
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We have measured the thermal conductivity of the spin S = 1 chain compound Y2BaNi0s. 
Analyzing the anisotropy of the thermal transport allows to identify a definite spin-mediated thermal 
conductivity k s along the chain direction. The calculated spin-related energy diffusion constant 
De(T) shows a broad peak around 120 K. Close to room temperature, De(T) approaches the 
theoretically predicted high-temperature value, while scattering of spin excitations by magnetic 
impurities seems to be the major limiting factor of n s at low temperature. 



PACS numbers: 66.70.+f, 75.40.Gb 75.47.-m 

A broad class of low-dimensional quantum spin sys- 
tems is expected to exhibit anomalous transport proper- 
ties due to the integrability of the corresponding model 
Hamiltonians [1, 2, 3]. A particularly well studied exam- 
ple is the one-dimensional Heisenberg antiferromagnetic 
(ID HAF) 5=1/2 XXZ model, for which several theo- 
retical works find finite Drude weights in the frequency 
spectra of the thermal k(uS) and the spin conductivity 
<j(u>) (see the recent review article [4] and the references 
therein). This, in turn, results in ballistic, dissipationlcss 
thermal and spin conductivities. These predictions are 
supported by experimental observations of enhanced en- 
ergy diffusion constants De and spin diffusion constants 
Ds in several physical realizations of the 5 = 1/2 ID 
HAF model [5, 6, 7, 8]. In contrast, in 5 = 1 HAF chains 
the first excited state is separated from the singlet ground 
state by an energy gap A and the system is not integrable 
[9]. This excludes infinite conductivities, but anomalous 
diffusion constants may still be expected [10]. Experi- 
mental results on the thermal transport in 5 = 1 chain 
compounds are scarce. The existence of spin-mediated 
energy transport in the 5 = 1 chain compound AgVP2S6 
has been reported [11], but the sample quality did not 
allow for unambiguously resolving whether the cause of 
the relatively low values of the spin-mediated thermal 
conductivity is of intrinsic origin or is caused by defects. 

In this Letter we report a study of the anisotropic ther- 
mal conductivity of large high-quality single crystals of 
Y2BaNiC>5, which represents a nearly ideal realization of 
the spin 5 = 1 HAF chain. The basic building blocks 
are chains of NiO@ octahedra with Ni 2+ ions (5 = 1) 
running along the a axis of the orthorhombic structure 
with lattice parameters a = 3.76 A, b = 5.76 A, and 
c = 11.33 A(space group Immm). The ID character of 
the magnetic properties is well established by susceptibil- 
ity and inelastic neutron scattering measurements yield- 
ing J w 250 — 280 K, A « 100 K, and \ J±/J\ < 10~ 4 for 
the in-chain superexchange coupling, the Haldane gap, 
and the ratio between out-of-chain and in-chain coupling, 
respectively [12, 13, 14, 15]. Our thermal conductivity 
data provide clear evidence for a sizeable spin-mediated 
energy transport along the chain direction, which is ab- 



sent along both perpendicular directions. In agreement 
with theoretical expectations, the derived energy diffu- 
sion constant shows a qualitatively different temperature 
dependence as those found for various S — 1/2 chain 
compounds. Our data also suggest that magnetic impu- 
rities provide the dominating scattering mechanism for 
magnons at low temperatures. 

The single crystals of Y2BaNi05 were grown by a self- 
flux method [16] . Out of one crystal, a bar-shaped sample 
was cut with dimensions 0.89 x 0.79 x 0.77 mm 3 along the 
principal axes a, b, and c, respectively. For the k(T) mea- 
surements we used a standard steady-state method [17]. 
The thermal gradient of AT ~ 0.1 K has been determined 
using a differential Chromel-Au+0.07%Fe thermocouple. 
Unaccounted heat losses due to the connecting leads and 
the sample surface radiation are estimated to be negligi- 
ble below about 200 K although at higher temperatures 
they can lead to somewhat overestimated values of k. For 
each orientation a few k(T) runs have been made with 
removing and re-attaching the thermal contacts to the 
sample. Different runs gave qualitatively similar temper- 
ature dependencies, but different absolute values varying 
up to ±10% for each orientation. This mainly arises from 
the uncertainty in determining the exact distance of the 
thermocouple junctions. Up to 14 T, wc did not find any 
magnetic-field dependence of the thermal conductivity. 

In Fig.l we present typical results for k 1 (T) (i = a, b, c) 
of Y 2 BaNi0 5 . Near 20 K, k 1 (T) for all three crystal di- 
rections exhibits a low-temperature peak, which is typical 
for insulators where thermal transport is of phononic ori- 
gin. The peak signifies a gradual change of the dominant 
phonon scattering mechanism from sample boundaries 
and extended defects at low temperatures to phonon- 
phonon scattering at high temperatures [18]. Another 
feature which is observed only for the a direction is 
a shoulder-like feature between about 80 and 140 K. 
Such an anisotropy has been observed in various low- 
dimensional spin systems such as 2D 5=1/2 materials 
[17, 19, 20, 21, 22, 23, 24], 5 = 1/2 ladder compounds 
[25, 26, 27], and 5=1/2 chain compounds [7, 28, 29, 30]. 
The conventional explanation is the existence of an ex- 
tra contribution to the thermal conductivity, beside the 
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FIG. 1: Thermal conductivity of Y2BaNi0s along (a axis) 
and perpendicular (b and c axis) to the spin chains. The solid 
(dashed) lines represent calculations of the phonon contribu- 
tions with (without) taking into account a minimum phonon 
mean free path. The inset shows the spin contribution to the 
thermal conductivity along the chain direction (see text). 



phonon conductivity (ft p h), coming from spin excitations 
(k s ), such that along the direction(s) of strong magnetic 
exchange = K p h + n s . Because of weak interactions 
perpendicular to the chains or planes, k s is negligible in 
these directions, that is, n 1 - ss K p h- However, one has to 
be careful in associating any high-temperature anomaly 
in kJI (T) with the spin transport, because similarly look- 
ing features can also appear in a purely phononic thermal 
conductivity as a result of scattering via spin-lattice in- 
teraction in certain temperature regions [31], as well as 
being caused by additional thermal transport via opti- 
cal phonon modes [32]. There are no obvious reasons, 
however, for these purely phononic mechanisms to show 
up only along the spin-chain direction. Thus, the most 
likely origin of the high-temperature feature in K a (T) of 
Y2BaNiC>5 is the spin contribution n s . 

In order to analyze k s (T) one has to subtract the 
phonon contribution K® h (T) from the total measured 
K a (T) along the a direction. Assuming that the same 
scattering mechanisms influence K 2 h along all directions, 



TABLE I: Parameters of the fits of k*(T) to Eq. (1). 



Parameter a axis b axis c axis 



a (i(r 17 s 1 ) 


4.0 


3.8 


2.4 


B (HT 18 s KT 1 ) 


9.1 


12.0 


14.9 


b 


6.7 


8.1 


6.1 


^min {^) 




18.3 


12.8 


Temperature region (K) 


5-40 


5-300 


5-300 



we first analyzed K b (T) and k c (T), assumed to be purely 
phononic, in order to establish the main phonon scat- 
tering mechanisms. Then we applied the same analysis 
to n a (T) at low temperatures where the spin contribu- 
tion k s (T) is expected to be negligibly small due to the 
spin gap, and extrapolated K^ h (T) to higher tempera- 
tures. For this analysis we employed the standard model 
of phonon thermal conductivity [18]: 



where x = fiuj/k^,T 1 u> is the frequency of a phonon, 
t(u,T) is the relaxation time, 6_d is the Debye temper- 
ature, v = 6/5 {ks/K) (67r 2 n)~ 1//3 is the mean sound ve- 
locity, and n is the number density of atoms. Due to the 
low k(T) ~ 1 W m -1 K _1 at high temperatures, which is 
comparable to n values of amorphous materials [33], we 
took into account a lower limit for the phonon mean free 
path £ m - m . The underlying idea is that the phonon mean 
free path cannot be smaller than roughly the size of the 
unit cell. Such an ^ m ; n is also observed in crystalline ma- 
terials at high temperatures [34]. Thus, the equation for 
r was put in the form t(u,T) = max{rs(w, T), £ m i n /v}. 
We now assume that all phonon scattering mechanisms 
act independently, such that t^ 1 = J^t^ 1 , where each 
term rf 1 corresponds to an individual scattering mecha- 
nism. Wc found that two scattering mechanisms, namely 
phonon scattering by extended two-dimensional defects 
and phonon-phonon scattering, in the form 

t s x = Alo 2 + Bu 2 Texp{-G D /bT), (2) 

are sufficient to approximate n b (T) and n c (T). We used 
On = 390 K estimated from the low-temperature specific 
heat [35]. The parameters A, B, b, and £ m i n were treated 
as freely adjustable parameters; the fit values are pre- 
sented in Table I. The fits of K ph (T) for i = b and c are 
shown in Fig. 1 (b,c) as solid lines. The dashed lines are 
obtained for the same values of A, B, and b, but setting 
^min = 0. This demonstrates that £ m i n does not influence 
Kp h (T) below about 100 K. Hence, only the parameters 
A, B, b can be adjusted by the fit of Kp h (T) in the tem- 
perature range T < 50 K. We do not attribute much 
significance to the absolute values of the fit parameters 
for two reasons: (i) the Debye model is a crude estimate 
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of the phonon spectrum of a real material, and (ii) the fit 
is only used to estimate the phonon background contri- 
bution, which has to be subtracted in order to separate 
the spin-mediated thermal conductivity. 

The spin thermal conductivity along the chains 
k s (T) = n a — Kp h is plotted in the inset of Fig. 1 (a). For 
Kp h (T) we used the solid line of Fig. 1 (a) which is calcu- 
lated for the fit parameters of Table I and ^ m in = 15.5 A, 
which is the average of the £ m ; n values for the b and c axes. 
The dotted lines are obtained for /tp h (T) if the individual 
values of £ m { n for b and c are used. One can see that the 
errors of n® h due to the uncertainty of £ m i n are negligi- 
ble below about 140 K and reach about ±10% at 300 K. 
Below about 70 K, k s is not well resolved because the dif- 
ference between the total measured k" and the estimated 
phonon contribution becomes of the order of relative 
experimental errors, therefore we will consider only the 
k s data above this temperature. The maximum of k s (T) 
amounts to about 5 W K _1 m -1 . That is about five times 
higher than that observed for the only other Haldane sys- 
tem AgVP 2 S 6 (J ~ 780 K;A- 240 K) for which thermal 
conductivity has been investigated so far [11]. From our 
k s (T) data, we calculated the spin-related energy diffu- 
sion constant D E (T) = K s (T)/(C s (T)a 2 ), where C S (T) 
is the specific heat of the spin chain and a the lattice con- 
stant along the chain. For C S (T), the numerical results 
for the ideal HAF 5 = 1 chain with J/k B = 280 K, cal- 
culated within the quantum nonlinear er-modcl (NLerM) 
[36], were used. 

The results of the present study together with the data 
for AgVP2S6 are shown in Fig. 2 (a). It is remarkable 
that not only the absolute values, but also their temper- 
ature dependencies are similar for both materials. An 
important observation is that, for both 5=1 materials, 
De(T) remains rather low in the whole investigated tem- 
perature region. This is in stark contrast to the S = 1/2 
chains where De(T) increases by about two orders of 
magnitude with decreasing temperature, as is also shown 
in Fig. 2 (a). 

The question is whether this is a generic difference, 
related to the integrability of the 5 = 1/2 and the non- 
integrability of the 5 = 1 chain model, or is it caused by 
extrinsic influences. To our knowledge, theoretical cal- 
culations for the energy diffusion in 5 = 1 HAF chains 
at T < J/ks are absent. The early analysis by Hu- 
bcr et al. gives the high-temperature limiting value for 
D E l ~ 2J/K [37], while a recent numerical study sug- 
gests a higher value of D E l « 5.6 J/h [10]. Indeed, at 
high temperatures De{T) of Y2BaNi05 is close to this 
limiting value, see Fig. 2 (a). The fact that De{T) ap- 
proaches D E l calculated for purely intrinsic interactions 
[10], suggests that intrinsic processes largely determine 
the energy diffusion in Y2BaNiC>5 at high temperatures. 

However, the experimental energy diffusion constant 
for 5=1 chains does not scale as T/J, although such 
a scaling should hold for intrinsic diffusion. This scaling 
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FIG. 2: (a) The energy diffusion constant De(T), calculated 
from the thermal conductivity data of Y2BaNiOr,. For com- 
parison, we also show De(T) of the 5 = 1 chain compound 
AgVP2S6 and of the 5 = 1/2 chain compounds BaCu2Si2C>7 
and SrCuC>2 (data from Refs. 7, 8, 11). The arrow corre- 
sponds to the high-temperature limit D 1 ^ calculated in [10]. 
(b) The normalized mean free path of the spin excitations cal- 
culated from the thermal conductivity data. The lower limit 
^mfiT 0nt °f the average chain segment length estimated from 
the susceptibility data is shown by the dashed line, (c) Illus- 
tration of the expected differences in defect scattering in the 
5 = 1/2 and 5 = 1 chains. The solid lines correspond to the 
mean free path for spin excitations when defect scattering is 
dominant. The dotted lines corresponds to a temperature- 
dependent scattering which is expected to become dominant 
at high temperatures, such as interaction between spin exci- 
tations or interactions with phonons. 



is also absent for the 5=1/2 case where predominantly 
external perturbations, such as spin-phonon interaction 
[39] , determine the behavior of spin-mediated heat trans- 
port. This suggests that, at least for a part of the excita- 
tion spectrum of the 5=1 chains, external perturbations 
are also important. Besides phonons, magnetic impuri- 
ties, that is ions with spin 5^1, are also expected to 
contribute to the scattering of spin excitations. If, in a 
certain temperature range, this type of scattering dom- 
inates, one may expect the magnon mean free path l s 
to be close to the average distance between the impuri- 
ties in the chains. Information about the concentration 
of magnetic impurities can be obtained from the low- 
temperature increase of the magnetic susceptibility, if the 
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type of impurities is known [40, 41]. For our sample, this 
is not known, but we can estimate the lower limit of the 
average chain segment length ^F™™* by assuming that 
all magnetic impurities have S = 1/2 and are situated 
in the chains. A corresponding fit of the susceptibility 
data (not shown) of our sample yields ^fiT™* ~ 60 A. 
In Fig. 2 (b) the mean free path of the spin excitations, 
estimated via £ s = k s /(C s v s ) where v s is the average ve- 
locity of spin excitations, is compared with ^fn"™*- In 
our calculations we used v s (T) = {ZK)~ X J e' k f(e,T)dk, 
where Z = / f(e,T)dk, f(e,T) = (exp(e/fc B T) - l)" 1 , 
and e(k) is the dispersion relation. For magnons in the 
5 = 1 chain we used the corresponding relations from the 
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A(T) 1 
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with A = 0.41J 



NLctM, e(k) = [V 2 {kaf 
and V = 2.49 J [38]. 

At low temperatures l s (T) approaches ^fn"™* su §" 
gesting that magnetic impurities are indeed one of the 
most important sources of magnon scattering. For sev- 
eral 5=1/2 chain compounds a saturation of £ S {T) at 
low temperatures is observed [8]. This is intuitively un- 
derstood as resulting from scattering by defects, in anal- 
ogy with the boundary scattering in conventional sys- 
tems. However, after reaching a peak around 120 K the 
derived £ S (T) of Y 2 BaNiOs decreases again with further 
decreasing temperature; a similar behavior was also ob- 
served for the 5 = 1 compound AgVP2Sg [11]. Yet, this 
behavior does not contradict the assumption that impu- 
rities are the dominant scatterers at low temperatures, 
since in gapped ID spin systems impurities are expected 
to give a constant mean free path only at high tempera- 
tures [42]. At T <C A, the scattering efficiency of point 
defects is energy-dependent and the low-energy spin ex- 
citations are scattered by impurities stronger than high- 
energy excitations [42]. This should lead to a mean free 
path decreasing with decreasing temperature for T < A 
in the 5 = 1 chains, in contrast to the gapless 5=1/2 
chains, as is illustrated in Fig. 2 (c). This expectation is 
in apparent agreement with our experimental results. 

In conclusion, our study of the anisotropic thermal con- 
ductivity of the 5 = 1 Haldane chain system Y2BaNiC>5 
permits an evaluation of the contribution of spin exci- 
tations to the heat transport along the chain direction. 
The corresponding energy diffusion constant is found to 
be close to the predicted high-temperature limiting value. 
With decreasing temperature, the spin-related energy 
diffusion constant of Y2BaNi05 increases much less than 
those of 5 = 1/2 chain compounds in agreement with 
theoretical expectations. The absolute values and the 
temperature dependence of the magnon mean free path 
suggest that magnetic impurities provide the dominating 
scattering mechanism for magnons at low temperatures. 

We acknowledge useful discussions with J. A. Mydosh. 
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